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1.  INTRODUCTION 

The  processing  of  digital  signals  is  rapidly  approaching  a  dominant 

»  * 

role  in  both  communication  and  radar  systems.  However,  the  physical 
origin  of  many  information-bearing  signals  is  intrinsically  analog  and 
continuous -time  in  nature.  Therefore,  an  effective  interface  between  the 
analog  and  digital  worlds  is  of  significant  importance  in  modern  signal¬ 
processing.  Very  often  the  quality  of  analog-to-digital  (A/D)  conversion 
is  the  critical  limiting  factor  in  over-all  system  performance. 

Quantization,  the  essential  mechanism  of  A/D  conversion,  begins  with 
the  availability  of  Input  values  ranging  from  -*  to  +•  .  The  quantizer 
replaces  each  of  these  input  values  with  a  prechosen  output  value.  The 
key  feature  is  that  each  output  value  is  one  of  a  finite  set  of  real 
numbers.  Hence,  a  symbol  from  a  finite  alphabet  can  be  used  to  represent 
and  identify  the  particular  output  values  that  occur.  Intuitively,  a 
different  n-bit  binary  word  can  then  be  associated  with  each  output  value 
if  the  set  of  output  values  contains  no  more  than  2n  members. 

The  moving  target  indicator  (MCI)  is  a  processor  required  in  most 
search  radars  to  separate  moving  targets  from  stationary  clutter. 
Nevertheless,  the  implementation  of  an  MTI  by  analog  means  has  been 
hindered  somewhat  because  of  the  instability  problems  incorporated  in 
an  MTI  system.  Quantization  of  an  MTI  system  with  digital  means  is  one 
of  the  possible  remedies  to  these  problems.  Because  the  characteristics 
of  a  quantizer  are  governed  by  the  choice  of  the  number  of  output  levels 
and  the  spacings  and  step  sizes  of  these  levels,  the  performance  of  an 
MTI  with  a  quantizer  will  depend  solely  on  the  above  parameters.  It  is 
the  effect  of  these  parameters  on  the  digital  Mil  that  will  be  discussed 


in  this  thesis. 


2.  A  BRIEF  INTRODUCTION  TO  MOVING -TARGET  INDICATION 
One  of  the  basic  requirements  for  present-day  tactical  radars  is 
desirable  performance  in  a  heavy  clutter  environment.  This  requirement 
often  leads  to  a  radar  system  approach  which  utilizes  a  moving  target 
indicator  (MTI)  to  extract  moving  targets  from  the  clutter  background. 

MTI  systems  operate  by  virtue  of  the  Doppler  spectral  differences 
between  returns  from  moving  targets  and  clutter.  For  instance,  if  the 
target  is  moving  relative  to  the  clutter  background,  the  MTI  can  be  applied 
to  filter  out  the  undesired  clutter  return  by  exploiting  the  differential 
Doppler  frequency  shift  produced  by  the  relative  target  to  clutter  radial 
motion  [1], 

2.1  Historical  Development 

In  early  radars,  most  MTI  systems  were  implemented  by  using  analog 
techniques.  In  the  1950' s,  MTIs  were  primarily  of  the  analog  single 
canceller  type  that  used  a  quartz  ultrasonic  delay  line  to  delay  the  radar 
return  by  a  time  corresponding  to  the  radar's  pulse  repetition  frequency 
(PRF) .  The  desired  MTI  effect  was  produced  by  subtracting  the  delayed 
return  from  the  present  return.  More  advanced  MTI  systems  employed  a 
double  canceller  that  was  formed  by  connecting  two  single  cancellers  in 
cascade.  This  configuration  not  only  produced  a  more  desirable  transfer 
function,  but  also  compensated  for  unbalances  in  each  of  the  individual 
cancellation  circuits.  Multiple  cancellers  are  discussed  and  analyzed 


by  Emerson  [21.  However,  it  has  been  shown  that  the  added  complexity  of 
the  multiple  cancellers  offers  little  advantage  in  terms  of  performance 
over  their  counterparts.  Thus  the  inability  to  control  the  delay  time  and 
the  system  gain  by  analog  means  have  made  it  hard  to  enhance  the  analog 
MTI  to  operate  toward  the  theoretical  limit. 

Also  during  the  1950's,  attention  had  been  focused  on  methods  to 
implement  better  MTI  operation  within  the  existing  hardware  constraints. 

It  was  found  that  applying  feedback  around  a  double  canceller  type  MTI 
allowed  excellent  control  of  the  shape  of  both  the  passband  and  rejection 
band  of  the  MTI  [3,4,5].  However,  when  this  technique  was  incorporated 
into  operational  radars,  a  significant  problem  was  observed:  any 
asynchronous  interference  present  caused  a  transient  that  circulated 
around  the  velocity  shaped  canceller  for  many  PRF  periods  resulting  in 
unacceptable  operation.  Hence  current  MTI  design  practice  tends  to  use 
nonrecursive  type  MTI  filters  wherever  possible  [6, 7,8, 9]  to  take 
advantage  of  their  superior  transient  response. 

More  recently,  because  of  the  availability  of  the  technology  of 
semiconductor  electronics  and  integrated  digital  devices,  microelectronic 
range-gated  filter  [2]  and  digital  MTIs  have  also  been  implemented.  The 
range-gated  filter  MTI,  as  the  name  implies,  range-gates  the  coherent 
MTI  signal  into  a  filter  that  separates  the  target  from  the  clutter. 
McAulay  [10]  has  shown  that  the  optimum  MTI  processor  is  achieved  when 
a  clutter  filter  is  cascaded  with  a  bank  of  contiguous  Doppler  filters 
matched  to  the  target  spectrum.  This  implementation  simultaneously 
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performs  Doppler  filtering  and  coherent  integration,  but  has  the 
disadvantage  that  a  large  number  of  filters  is  required.  On  the  other 
hand,  the  digital  MTI  is  basically  a  digital  implementation  of  the 
analog  MTI.  In  fact,  digital  techniques  are  not  only  practical  but  they 
are  essential  to  meet  the  present  day  tactical  requirements  of  light 
weight,  small  volume,  high  availability,  and  near-optimum  system  per¬ 
formance  at  a  reasonable  cost.  Hence,  the  relative  high  cost  of  range¬ 
gated  filter  MTIs  and  the  development  of  low  cost  digital  components 
results  in  a  general  preference  for  digital  MTIs. 

The  MTI  problem  for  surface-based  radars  is  less  complex  than  that 
for  airborne  radars  [2,11].  The  complications  in  airborne  radars  arise 
primarily  because  of  the  clutter  spectral  spreading  and  translation 
effects  caused  by  platform  motion  in  addition  to  the  aircraf t-to-surface 
geometry.  There  are  two  general  approaches  for  solving  the  airborne  MTI 
(AMTI)  clutter  problem.  The  first  is  generally  referred  to  as  a  pulse 
Doppler  radar.  Pulse  Doppler  radars  are  used  for  airborne  interceptor 
(AI)  radars  where  "look  down"  in  heavy  clutter  is  required.  One 
application  of  this  technique  is  to  the  AWACS  EC-121  airborne  radar  [12]. 
The  other  approach  for  solving  the  AMTI  radar  clutter  problem  is  to 
remove  the  effects  of  aircraft  motion  by  antenna  design  and  signal 
processing.  Once  the  platform  motion  effects  are  removed,  the  same 
signal  processing  techniques  as  used  in  the  surface  MTIs  are  applicable. 
Rhys  and  Andrews  [13]  has  a  detailed  description  of  this  type  of  system. 
Application  of  this  technique  is  to  the  E2-C  airborne  early  warning 
(AEW)  radar. 


2.2  Clutter  Modeling 

Before  we  analyze  the  performance  of  the  MTI  radars,  the  characteristic 
of  the  clutter  environments  should  be  discussed.  Since  an  MTI  system 
might  perform  better  in  one  clutter  environment  than  in  another,  a  clutter 
model  assumption  is  implicit  in  the  design  or  analysis  of  the  system. 

Ideally,  the  model  used  should  accurately  reflect  the  clutter  process 
in  a  realistic  manner.  Practically,  the  clutter  data  base  is  rarely 
adequate  to  achieve  this  objective,  and  the  art  of  modeling  the  clutter 
involves  a  compromise  between  accurate  representation  and  analytic  con¬ 
venience.  The  following  is  a  brief  overview  of  the  various  clutter 
models  used  commonly  in  the  design  or  analysis  of  the  MTI: 

(i)  The  Gaussian  (Rayleigh  Envelope)  Clutter  Model 

The  Gaussian  or  Rayleigh  envelope  clutter  model  is  applicable 
to  distributed  clutter  that  is  returned  from  a  spatially  continuous 
distribution  of  scatterers  with  no  subset  of  scatterers  predominating. 

This  type  of  clutter  is  usually  associated  with  weather  clutter,  chaff, 
sea  clutter  observed  with  a  low  resolution  radar  (pulse  width  >  0.5  sec) 
or  with  a  high-resolution  radar  at  high  grazing  angles  (cp  >  5°)  and  land 
clutter  observed  for  high  grazing  angles  (cp  >  5°)  over  undeveloped 
terrain. 

The  random  clutter  process  is  represented  at  the  clutter  carrier 
frequency  (f£)  as 


c  =»  x„  cos  m  t  -  y  sin  cu  t 
t  t  c  J  t  c 


(2.1) 


where  x  and  are  zero  mean,  identically  distributed,  low  pass, 
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independent  normal  processes  with  variance  a  .  The  clutter  process  c^. 
is  strict-sense  stationary  and  hence  its  characteristics  are  independent 
of  the  time  origin. 

The  power  spectral  density  associated  with  the  randomly  fluctuating 
portion  of  the  clutter  cross-section  is  usually  assumed  to  be  Gaussian 
shaped  in  that  it  is  the  result  of  several  independent  effects.  The 
power  spectral  density  of  the  power  envelope  is  given  by 

p 

Sn(f)  -  — =  exp[-f2/2a  2]  (2.2) 

P  /  J  f 

»/2TTaf 

where  is  the  clutter  power  and  a ^  is  the  standard  deviation  of  the 
spectrum. 

The  autocorrelation  function  associated  with  the  power  envelope  is 
given  by 

Rp(T)  -  Pc  exp[-T2/2CTT2]  (2.3) 

where  a  »  1/2ttu,. 

T  f 

The  normalized  autocorrelation  function  is  then 

Pc(T)  =  Rp(T)/Pc  =  exp[-T2/2aT2]  (2.4) 

(ii)  The  Rician  Clutter  Model 

The  Rician  clutter  model  is  similar  to  the  Gaussian  model 
except  that  a  predominant  steady  scatterer  (S)  is  added  to  the 
distributed  clutter. 

The  random  clutter  process  is  represented  by 

C  =*  (S-Hx„)  cos  cd  t  -  y  sin  U)  t  . 
t  t  c  t  c 
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The  autocorrelation  function  for  the  in-phase  process  and  Gaussian 
shaped  distributed  clutter  is 

R(t)  *  P0[m2  +  exp(-T2/2crT2)  ] 

2  2 

where  P  is  the  distributed  power  and  m  *  S  /P  . 
o  o 

The  resulting  power  spectral  density  is 

S(f)  =  m2PQ6(f)  +  PQexp(-f2/2CTf2)  /Jlua* . 

(iii)  The  Log-Normal  Clutter  Model 

The  log-normal  clutter  model  is  associated  with  high  resolution 
(pulse  width  <  .5  sec)  sea-clutter  data,  where  the  sea  clutter  is  observed 
at  grazing  angles  (9  <  5°);  and  ground  clutter  observed  at  low  grazing 
angles . 

The  random  clutter  process  is  represented  by 

/  2  2 

ct  =  ^Xt  +  yt  cos  (V  +  V 

where  cpc  is  a  uniformly  distributed  random  variable. 

The  log  normal  (2  )  variable  is  related  to  the  Gaussian  variable 
(xfc)  by  a  nonlinear  function,  A  exp(*)>  where  A  is  a  constant  and  in 
which 

R  (0) 

E(Zt)  =  A  exp  [-^2 — 1 
Rz(t)  -  A2exp[Rx(0)  +  Rx(T)] 

Cz(t)  =  A2exp[Rx(0)][exp[Rx(0)]-l3 

where  &X(T):  autocorrelation  function  (ACF)  of  the  Gaussian  variable 
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R^Cr):  ACF  of  the  log-normal  variable 
C  (t):  covariance  function  of  the  log-normal  variable. 

L  » 

(i v)  The  Weibull  Clutter  Model 

The  Weibull  clutter  model  is  used  to  model  both  land  and  sea 
clutter  and  offers  the  potential  to  accurately  represent  the  real  clutter 
distribution  over  a  much  wider  range  of  conditions  than  either  the  log¬ 
normal  or  Gaussian  model.  Szajnowski  [14]  has  the  derivation  of  the 
correlation  properties  of  the  Weibull  clutter  model.  He  showed  that  the 
correlation  function  of  the  Weibull  clutter  process  is  related  to  that  of 
the  Gaussian  clutter  process.  Moreover,  the  Weibull  (Wfc)  variable, 
which  is  governed  by  two  parameters  p  and  q,  is  related  to  two  independent 
Gaussian  variables  x^  and  y  such  that 

Wt  *  (xt2  +  yt2)1,/p 

where  xt  and  y^  are  zero-mean  identically  distributed  Gaussian  variables 
2 

with  variance  ct  ,  p  is  the  sample  parameter  and  q  is  the  scale  parameter 

2  2  1/p  2 

related  to  ct  by  q  =>  (2ct  )  .  The  mean  e£w  }  and  the  variance  ct  of 

c  w 

the  Weibull  variable  are 

E{Wt3  -  (2CT2)1/p  T(1  +  1/p) 
ctw2  *  (2ct2)2/p£(1  +  2/p)  -  r2(l  +  1/p)} 
where  T(*)  is  a  garrena  function. 
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Of  the  above,  the  simplest  clutter  model  is  the  Gaussian  model 
whereby  the  random  clutter  process  is  represented  as  a  Gaussian  process 
which  is  completely  specified  by  its  mean  and  correlation  function.  At 
present,  analytical  techniques  are  not  available  to  determine  the 
performance  of  MTI  systems  in  the  presence  of  non-Gaussian  clutter.  The 
principle  method  of  attacking  this  type  of  problem  is  through  a  Monte 
Carlo  computer  simulation  in  which  the  actual  MTI  system  is  exercised  by 
applying  statistical  samples  of  the  postulated  clutter  and  target 
processes.  For  this  reason,  the  Gaussian  model  has  been  used  exclusively 
in  the  MTI  analysis.  Henceforth,  we  shall  choose  the  Gaussian  model 
instead  of  the  others  as  the  clutter  model  throughout  the  remainder  of 
the  thesis . 

2.3  Mil  System  Performance 

The  ultimate  measure  of  the  MTI  system  performance,  i.e.,  the 
capability  of  the  MTI  to  separate  moving  targets  submerged  in  the 
clutter,  does  not  depend  on  power  considerations  alone.  The  following 
performance  criteria  are  seen  most  often: 

(i)  Clutter  Attenuation  (CA) 

This  criterion  is  defined  as  the  ratio  of  the  power  in  the 
clutter  before  Doppler  filtering  to  the  power  in  the  clutter  after 
Doppler  filtering,  with  the  output  clutter  power  normalized  to  unity 
gain  for  uncanceled  signals;  i.e.. 


CA 


Sc(m)dm 

Sc(<u)lH(u))|2diu 
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where  Sc(u))  Is  the  clutter  power  spectral  density  and  H(<u)  is  the  MTI 
filter  transfer  function. 

(ii)  Subclutter  Visibility  (SCV) 

Subclutter  visibility  is  defined  as  the  ratio  at  the  input  of 
the  Doppler  filter  of  clutter  power  to  moving  target  power  to  produce 
equal  Doppler  filter  outputs.  Specifically, 

C  S_(u)-<ud)  (H(u))  |2duj 
SC V  -  -  ;  pc  -  pt 

Sc(u))|H(u>)|2du> 

where  ST(cu-(ud)  is  the  target  power  spectral  density.  The  SCV  equation  is 
usually  evaluated  at  <u^  equal  to  an  optimum  speed  target;  i.e.,  a  target 
whose  Doppler  frequency  is  one-half  the  radar's  pulse  repetition  frequency. 
However,  it  is  more  fundamental  to  evaluate  the  SCV  at  the  target  speeds 
of  interest. 

(iii)  Improvement  Factor  (or  Reference  Gain)  (I) 

This  criterion  is  found  by  comparing  the  input  signal-to-clutter 
ratio  to  the  Mil  output  signal-to-clutter  ratio  for  targets  distributed 
uniformly  over  the  radial  velocity  spectrum. 

(Wo 

I  »  .  ;  average  over  all  possible  target 

's'  c'i  Doppler  frequencies 

(iv)  Mil  Gain 

The  gain  of  an  Mil  is  defined  as  the  output  of  the  Doppler 
filter  when  excited  with  a  process  whose  spectral  density  corresponding 
to  all  velocities  of  interest  is  uniformly  distributed. 


J“max  Jh  (u))|2du> 

MTI  gain  -  — - - - 

°  OJ 

max 

2 

where  u>mav  is  the  maximum  Doppler  frequency  of  interest,  and  (h^Cuu)  J 
is  the  normalized  power  transfer  function. 

Clutter  attenuation  and  subclutter  visibility  are  equivalent 
criteria  when  target  and  clutter  spectral  shapes  are  similar  and  the 
response  of  the  Doppler  filter  is  unity  for  the  targets  under  considera¬ 
tion.  For  single  and  double  cancellers,  this  occurs  for  optimum-speed 
targets.  The  improvement  factor  provides  a  criterion  for  computing  the 
average  improvement  in  the  signal- to-clutter  ratio  introduced  by  the 
canceller;  whereas.  Mil  gain  provides  a  criterion  for  comparing  target 
responses  of  various  MTI  filters. 

In  general,  the  criterion  that  best  represents  the  desired  objective 
of  the  system  should  be  selected.  For  example,  consider  an  MTI 
surveillance  requirement  where  the  moving  target  is  to  be  located  by 
examining  a  display.  The  MTI  system  must  reduce  the  ground  clutter  to 
a  level  consistent  with  the  radar's  noise  level  to  keep  the  false 
alarms  at  a  tolerable  level,  while  allowing  the  detection  of  weak 
targets.  For  this  type  of  situation,  the  clutter  attenuation  and  the 
MTI  gain  are  appropriate  performance  measures;  the  first  relating  to 
the  false-alarm  performance,  the  second  to  the  ability  to  detect  a 
variety  of  targets.  As  a  second  example,  consider  a  combat  surveillance 
radar  whose  function  is  to  detect  personnel  whose  velocities  are  1  to  3 


12 


knots.  The  performance  criterion  of  interest  here  is  maximum 
subclutter  visibility  evaluated  at  1  to  3  knots  velocities.  For 
our  discussion,  the  quantization  effects  on  the  digital  MT1  are  to 
be  analyzed,  the  desirable  performance  measure  in  this  case  is  the 
clutter  attenuation  since  it  is  most  sensitive  to  the  errors  introduced 
by  the  quantizer  [ 15 ] • 
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3.  PERFORMANCE  EXPRESSIONS  FOR  DIGITAL  MTI 
3.1  An  Analytical  Model  for  Digital  MTI 

The  performance  of  analog  MTI  systems  has  been  limited  by  problems 
in  controlling  the  delay  times  and  gains  of  such  systems.  However,  these 
problems  have  been  overcome  by  a  digital  MTI  system  which  is  basically  a 
system  with  a  quantizer  being  inserted  at  the  input  of  the  analog  MTI 
system.  In  the  following  discussion,  the  previous  works  related  to  the 
quantization  effects  on  the  digital  MTI  will  be  reviewed;  however,  the 
primary  emphasis  will  be  on  the  single  and  double  canceller  MTIs. 

Mathematically,  the  output  y(t)  of  a  single-canceller  MTI  is  related 
to  the  input  voltage  x(t)  by 

y(t)  -  x(t)  -  x(t+D)  (3.1) 

where  D  is  the  system  delay;  and  similarly  for  the  double  canceller  MTI, 
the  output  is 

y (t)  -  x(t)  -  2x(t-H>)  +  x(t+2D)  (3.2) 

Eqs.  (3.1)  and  (3.2)  are  for  linear  MTIs.  If  a  quantizer  is  inserted 
before  MTI  delay,  then  outputs  of  the  single  and  double  canceller  digital 
MTIs  are 

yQ(t)  -  Xq(t)  -  Xp(t+D)  (3.3) 

and 

yQ(t)  -  xQ(t)  -  2xQ(t+D)  +  Xq  (t+2D) 
respectively,  where  x^(t)  is  a  quantized  version  of  the  input  x(t). 


(3.4) 


The  block  diagrams  for  a  single  canceller  MTI,  a  double  canceller 
MTI  and  a  digital  MTI  are  shown  in  Figs.  3.1,  3.2,  and  3.3,  respectively. 


3.2  Distortion  Effects  Due  to  Digitizing  MTI 

Before  deriving  the  equations  needed  for  the  analysis  of  the  Mil 
performance,  we  will  examine  the  distortions  introduced  in  the  presence 
of  a  quantizer. 

In  digital  MTI,  the  quantizer  is  the  critical  element  of  the  design, 
as  the  ultrasonic  delay  line  is  to  the  analog  MTI.  Because  of  the 
presence  of  the  quantizer,  digitized  input  generates  basically  two  types 
of  distortion. 

The  first  is  the  quantization  noise  that  must  be  considered  in  addition 
to  the  usual  noise  encountered  in  the  analog  MTI.  Brennan  and  Reed  [16] 
first  defined  the  output  ratio,  Rq,  as  the  performance  criterion  of  an 
MTI  system,  where 

|ECyQ(t) | signal  present}  -  E[yq(t) (signal  absent}|2 
Ro  var{yq(t) | signal  absent] 

The  numerator  is  the  square  of  the  incremental  digitized  output  voltage 

in  the  presence  of  a  signal,  and  the  denominator  is  the  variance  of  the 

digitized  output  voltage  in  the  presence  of  the  clutter  plus  quantization 

noise.  In  general,  the  expectation  of  y^(t)  is  zero  in  the  absence  of 

signal  [16].  Thus  Rq  is  simplified  to  the  form  in  which 

2 

(E{y  (t)( signal  present] ( 

^  ^  —  — 

°  E { | y^ ( t ) | 2| signal  absent} 


Delay 

D 


Figure  3.1.  Single  canceller  MTI 


1 

Delay 

D 

1 

Figure  3-2.  Double  canceller  MTI 


Figure  3-3.  Digital  MTI. 
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Brennan  and  Reed  also  identified  a  quieting  effect  due  to  correlation 
between  quantization  errors  on  successive  pulses.  For  the  case  of  interest, 
assume  a  uniform  quantizer  is  utilized,  the  output  ratio  for  a  single 
canceller  Mil  is  then  (recall  that  we  are  assuming  Gaussian  statistics 
throughout) 


<Vsi>‘ 


2ax2[ !-px(D) ]  +  ~  ( 1— Q) 


(3.5) 


2 

where  •  E[x^(t+iD) j signal  present];  is  the  variance  of  the  input 
process  x(t);  p^(D)  is  the  normalized  correlation  function  of  the  input 
process  x(t) ;  is  the  quantization  interval  and  Q  is  a  quieting  factor 
which  is  a  measure  of  quieting  due  to  correlation  between  quantization 
errors  on  two  successive  pulses;  and  mathematically  Q  is  given  by 


Q  =  ~f  2  exp[-4TT2k2(o  /A)2[l  -  P  (D)  ] } 
TT  k=l  k  XX 

For  a  double  canceller  MTI ,  the  output  ratio  is 
„ _ (Sp  -  2St  -  S2)2 


(3.6) 


°  2  A2 

2ax  [3_4px(D)  +  PxC2D)  ]  +  ^(1-Q) 


(3.7) 


In  this  case,  the  quieting  factor  is 


Q  -  -|  2  -j  exp{-4TT2k2((a  /A)2[l-p  (D)]} 
tt  k-1  k  xx 


+  ~§  2  “4  exp[-4TT2k2(cr  /A)2[l-p  (2D)  ]} 

tt  k-1  k  xx 


(3.8) 
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Each  of  the  denominators  of  Eqs .  (3.5)  and  (3.7)  contains  two  terms.  The 
first  is  the  residue  clutter  power  in  the  absence  of  quantization  noise, 

ft 

and  the  second  is  the  noise  power  due  to  quantization.  Without  the 

quieting  effect,  as  for  the  single  canceller  MTI,  the  second  term  simplifies 
2 

to  A  /6,  which  can  be  obtained  otherwise  by  treating  the  quantization  errors 
in  the  Gaussian  process  x(t)  as  independent  and  uniformly  distributed 
between  -A/2  and  A/2. 

The  second  is  the  clipping  error  which  occurs  when  the  input  voltage 
is  greater  than  the  maximum  quantizer  level.  The  dynamic  range  of  the  input 
is  seriously  restricted  in  the  presence  of  the  quantizer,  a  form  of  limiter, 
which  in  turn  regulates  the  clutter  residue  at  the  digital  MTI  output. 
Unfortunately,  the  quantizer  causes  the  spectrum  of  the  clutter  to  spread 
so  that  the  MTI  performance  in  the  presence  of  limiting  falls  short  of  the 
performance  predicted  from  linear  theory.  Ward  and  Shrader  [17]  computed 
the  MTI  improvement  factor  from  the  autocorrelation  function  for  limited 
clutter  (i.e.,  clutter  after  bandpass  limiting).  They  chose  a  clutter 
model  with  a  Gaussian  autocorrelation  function,  and  a  limiter  with  an  error 
function  characteristic  because  a  simple  expression  derived  by  Baum  [18] 
relates  the  input  and  limited  output  autocorrelation  functions  when  the 
input  statistics  are  Gaussian.  The  results  they  found  are  verified  by 
time-domain  Monte  Carlo  simulation.  They  concluded  that  limiting  must 
be  included  in  any  analysis  of  the  MTI  performance  in  order  to  achieve 
meaningful  results.  The  presence  of  limiting  may  degrade  the  MTI 
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improvement  factor  by  a  large  amount.  For  example,  in  the  analysis  of  the 
double  canceller  MTI,  a  degradation  as  large  as  20  dB  may  occur  [17].  In 
the  presence  of  another  clutter  model,  in  which  the  input  statistic  is 
neither  Gaussian  nor  quasi -Gaussian,  Grasso  [19]  has  shown  that  limiting 
has  the  similar  impact  on  the  degradation  of  the  MTI  performance.  He 
summarized  that  the  improvement  factor  remains  practically  unchanged  in 
the  case  of  a  single  canceller  MTI,  but  that  the  improvement  factor  of 
the  double  canceller  is  seriously  reduced. 

Tong  [15]  analyzed  the  performance  of  the  MTI  with  uniform  quantizers. 
In  his  paper,  he  computed  the  clutter  attenuation  of  the  MTI  by  making  use 
of  an  expression  derived  by  Velichkin  [20]  which  utilizes  the  quantized 
normalized  correlation  coefficient.  This  expression  takes  into  account 
the  combined  effects  of  both  types  of  errors  discussed  above.  Since  we 
are  going  to  rely  heavily  on  this  expression  for  obtaining  the  MTI  clutter 
attenuation,  we  shall  derive  it  here  for  completeness. 

3.3  Output  Correlation  and  Clutter  Attenuation 

Given  a  zero-mean  Gaussian  stationary  random  process  x(t),  it  is 
level-quantized  by  a  symmetrical  quantizer  with  a  step  characteristic 


Q[x(t)]  which  is  given  by 


M-l  x(t) 

Q[x(t)]  »  +  2  J  Lyp  Cn-ak) dT) 

k=l  -<*> 


(3.9) 


where  5  (H-a^)  is  the  delta  function,  M  is  the  total  number  of  quantization 
levels  of  the  step  characteristic;  w^  is  the  first  output  level; 

■  w^^-w^  ts  the  k  step  size  of  the  quantizer;  a^  is  the  k  threshold 
spacing  of  the  quantizer.  The  typical  quantizers  are  shown  in  Figs. 

3.4  and  3.5. 


The  autocorrelation  function,  Rq(t)>  of  the  quantized  Gaussian  process 
is  the  joint  moment  of  Q[x(t-K)]  and  Q[x(t)]  and  is  defined  as 

00  00 

Rq(t)  =  E[Q[x(t4r  )  ]Q[x(t)  ]}  =  f  J  Q(u)Q(v)f  (u,v;t  )dudv  (3.10) 


where  u  =  x(t-K);  v  =  x(t);  f(u,v;r)  is  the  second  order  Gaussian 
probability  density  of  the  process  x(t). 

The  second  order  Gaussian  probability  density  can  be  written  in  terms 
of  the  Hermite  polynomials  and  the  corresponding  first  order  probability 
densities  [21] 


“  P°(  T) 

f(u,v,T)  =>  f(u)f(v)  £  Hn(^)HnC^) 

n=0  '  x  x 


(3.11) 


where  f(x)  is  the  first  order  Gaussian  probability  density;  PX(T)  is  the 
normalized  correlation  coefficient  of  x(t) ;  and  Hn(x)  is  the  Hermite 
polynomials  of  degree  n  defined  by  [22] 


2  ,nr  -x  /2\  2,„ 

Hn(x)  =  (-l)neX  2  d-~ =  (-l)11^  eX  2  $(n+1)i 

dx 


(3.12) 


where 


*  (t)  *  J  —  e'X  /2dx 
-<*>  V*2tt 


(3.13) 


Substituting  (3.12)  into  (3.11),  f(u,v,T)  becomes 

03  n(  \ 

1  VT  ,  (n+1)  ,  u..(n+l),  v. 

f  (u ,  v ,  T )  =  — 2  ^  ~ <5“) 


o  n=0 
x 


(3.14) 


Putting  (3.14)  into  (3.10),  Rq(z)  is  then  written  as 
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1  *  Pxn<T)  f 

RQ^T>  "  “I 


a  *  n*0  “*  L 


f  Q(u)f(n+1)(^)du'2 


(3.15) 


Since  Q(u)  is  odd  symmetric,  and  (~)  is  even  symmetric  for  n  even 

x 

the  sum  of  the  terms,  where  n  is  even,  equals  to  zero.  With  this 
mocification,  Eq.  (3.15)  becomes 

PXV)  r 

nl 


VT>--7  =  ^!i"Q(u)f(n+1)(^)du:2 


2  . 
a  n=l 

X  n  odd 


(3.16) 


Using  the  technique  of  integration  by  parts 


J  Q(u)$(n+1)(^)du  =  -  CTx  J*  Q(1)(u)$(n)(^i)du 


(3.17) 


Differentiating  (3.9)  with  respect  to  u,  we  get 


Q(1)(u)  =  2  Ak§(u-ak) 
k=l 


M-l 


(3.18) 


Substituting  (3.18)  into  (3.17),  we  obtain 

M-l 

2 

k*l 


f  Q(u)»<“+»(Ji)du  »  -  V 


(3.19) 


With  the  help  of  Eq.  (3.19),  R  (t)  has  the  form 


00  fM-1  /  \  *1,  ^2  Px 

VT)  -  z.  i,£,  v  <r>:  “t; 


n*l  t  k=l 
n  odd 


(3.20) 


;(n) 


Since  ?  (x)  is  related  to  H^_^(x)  by  Eq.  (3.12),  Eq.  (3.20)  becomes 


»  r  M-l  ,  ,  -2  p  "<T) 

VT)  ‘  \  J Vxp<-\  n\  VV ,  asn— 

n=  i  _  k=  i 
n  odd 


(3.21 


The  variance,  ,  of  the  quantization  process  is  determined  by  equating  r 
to  zero  in  Eq.  (3.21).  Then 


aQ2  -  V0)  .  £  [  V  Vxp<-ak2/2cx2)H[i.1(.k/0x),y 

n  odd 


(3.22 


The  normalized  correlation  coefficient,  Pq(t)  of  the  quantized  process  is 


VT)  Ro<T>  1  “r11"1  2  2 

P0(T)  *  R(oT  *  2“  =  - 2  S  2  Vxp<_ak  /2CTx  )Hn-l(ak/ax) 

y  oQ  2m  n=l  L  k=l  x  K  x  n  1  k  x 

"  ^  n  odd 


P  n(T) 
Kx 


(3.23 


Since  Pq(t)  and  ct^  are  infinite  series,  they  can  only  be  evaluated 
numerically.  A  stopping  rule  for  the  evaluations  of  them  is  provided  in 
Tong ' s  paper  [ 15 ] . 

The  clutter  attenuation,  a  performance  criterion  of  an  MTI  system, 

relates  the  input  power  and  the  output  power  of  the  MTI.  Mathematically, 

for  the  Gaussian  clutter  model  in  which  the  input  clutter  x(t)  is  a  zero 

2  2 

mean  stationary  Gaussian  random  process,  the  input  power  is  *  Etx^(t)} 

2 

The  output  power  ay  ,  of  a  single  canceller  MTI  is 


<7y2  “  E[y2(t)}  =*  E[[x(t)-x(t+D)  ]2}  =  2ct..2[1-p..(D)] 


(3.24 
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and  the  output  power  of  a  double  canceller  MTI  is 


a  2  -  EC[x(t)-2x(t+D)+x(t+2D)]2}  -  60^1-  J  PX(D)+  J  Px<2D>] 


(3.25) 


The  clutter  attenuation,  CA,  for  these  MTI's  are 


Ld>  =  a  _ i _ 

a2 

y 


,(2)  _ 


6[1  -  |  Px(D)  +  j  Px(2D>] 


respectively.  Similarly,  the  CAs  for  the  quantized  MTI  system  are 


2[ 1-Pq(D) ] 


k  (2)  _ 1 _ 

Q  6[l-|»Q(I»+ipQ<2D)J 


(3.26) 


(3.27) 


(3.28) 


(3.29) 


Notice  that  for  the  double  canceller  MTI,  CAV  7  is  a  function  of  PX(D) 
and  px(2D) .  However,  Px(°)  and  px(2D)  are  directly  related  to  one 
another  for  the  Gaussian  assumption.  Recall  the  correlation  function  for 


the  Gaussian  clutter  model  is  given  by  Eq.  (2.4)  and  is  repeated  here 


Px(T)  -  exp(-T2/2aT2)  . 


(3.30) 


Similarly  Px(2D)  can  be  shown  as 

p  (2D)  -  exp(-(2D)2/2a  2)  »  [exp(-D2/2a  2)]4  -  p  4(D).  (3.31) 

X  •  T  X 
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Hence,  the  clutter  attenuation  for  both  linear  MTIs  are  described  solely 
by  the  linear  correlation  coefficient,  PX(D),  while  the  clutter  attenua¬ 
tion  for  both  quantized  MTIs  are  governed  by  the  two  corresponding 
quantized  correlation  coefficients  p^(D)  and  p^(2D),  which  are  both  in 

terms  of  the  step  size  the  threshold  spacing  a^,  the  variance  of  the 
2 

input  process  ,  the  total  number  of  the  quantization  levels  M,  and  the 
linear  correlation  coefficient  p^(D). 

In  order  to  compare  the  performance  of  the  linear  MTI  and  the 
quantized  MTI,  the  loss  factor,  y>  is  defined  as  the  loss  introduced  due 
to  quantization,  in  which  for  the  single  and  double  canceller  MTIs  are 


and 


CA 


(1) 


CA, 


(1) 


Y  = 


CA 


(2) 


CA, 


(2) 


(3.32) 


(3.33) 


respec tively . 

In  the  following  sections  of  this  thesis  Eqs.  (3.23),  (3.3) ),  (3.32), 
and  (3.33)  will  be  used  extensively  for  analyzing  the  performance  of  the 
single  and  double  canceller  MTIs  with  uniform  quantizers,  logarithmic 
quantizers  and  p,-law  quantizers  respectively. 

3.4  Uniform  Quantization  in  Digital  MTI 

With  the  above  tools,  Tong  [15]  has  analyzed  the  performance  of  the 


single  and  double  canceller  MTIs  with  uniform  quantizer  as  a  function  ot 
the  numoer  of  quantizer  levels  and  the  spacing  of  these  levels.  His 
analysis  considers  only  the  even  quantizer,  i.e.,  M  is  even.  A  uniform 


quantizer,  as  its  name  implies,  has  a  uniform  threshold  spacing  and 
uniform  step  size.  Without  loss  of  generality,  unit  step  size  is  assumed. 
For  simplicity,  a  normalized  spacing  between  thresholds  of  an  M-level 


uniform  quantizer  is  denoted  as  a  *  a^/c^.  ^he  uniform  threshold  spacings 
of  the  uniform  quantizer  are  then  located  at  tea  where  k  is  an  integer 
which  goes  from  -(M/2-1)  to  (M/2-1)  including  k  *  0.  With  the  above 
modifications,  Pq(t),  for  the  uniform  quantizer,  can  be  shown  to  be 


r  M/2-1 


Po(T>  =  2 

"  to  n-1  L  k*0 
Q  n  odd 


2  &kexp(-(ka)  /2)Hn_1(ka) 


2  Pxn(T) 


(3.34) 


where  A^  *  \  for  k  *  0,  and  A^  ■  1  for  k  f  0. 

In  order  to  test  our  computer  program,  we  verified  Tong's  results  by 
using  Eq.  (3.34)  for  uniform  quantizer  with  even  numbers  of  levels. 

Examples  of  the  results  are  given  in  Figs.  3. 6-3. 9  which  depict  the  loss 
of  the  single  and  double  canceller  MTIs  with  8  and  32  level  uniform 
quantizers.  On  the  basis  of  similar  results,  Tong  concluded  that  the 
best  uniform  quantizer  threshold  spacing  is  a  function  of  the  input  clutter 
power  and  the  clutter  correlation  property.  The  best  a  is  found  to  be 
strongly  dependent  on  the  value  of  p^,  especially  in  the  case  when  the 
number  of  quantizer  levels  is  small  (see  figs.  3.6  and  3.8).  There  is 
always  an  optimum  a  for  each  input  correlation  function  p^(D).  A 
nonoptimum  spacing  can  result  in  a  loss  significantly  greater  than  the 
minimum.  Therefore,  if  clutter  input  is  expected  to  change  widely, 
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optimum  threshold  is  desirable  to  adapt  to  local  clutter  conditions. 


Furthermore,  as  the  number  of  output  levels  is  increased,  the  minimum 
loss  is  reduced  considerably  and  the  optimum  a  setting  is  less  critically 
dependent  on  the  input  correlation  coefficient.  However,  if  the 
threshold  spacing  is  either  toe  small  or  too  big,  as  Tong  noted,  the 


quantizer  acts  as  a  hard  limiter  and  the  loss  is  maximum, 
situation,  is  a  function  of  p  and  is  given  by  [23] 


in  this 


„  2  ■  -1 
PQ  =  TT  310  PX 


(3.35) 


For  instance,  when  p  is  0.999,  p  approaches  0.9715,  and  '{  approaches 

x  y  * 

14.5  dB  and  39.8  dB  for  the  single  and  double  canceller  MTIs  respectively. 

By  using  techniques  similar  to  those  of  Tong,  the  performance  of  the 
MTI  with  logarithmic  and  n-law  quantizers  will  be  discussed  in  the 
following  section. 
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4.  DIGITAL  MTI  PERFORMANCE  FOR  M-LAW  AND  LOGARITHMIC  QUANTIZERS 


In  this  section,  the  effect  of  logarithmic  and  p.-law  quantization  on 
the  performance  of  the  single  and  double  canceller  MTI  will  be  analyzed. 
Only  odd-symmetric  quantizers  with  number  of  output  levels  equal  integral 
powers  of  2  will  be  considered.  Without  loss  of  generality,  unity  step 
size  and  unit  variance  of  input  clutter  process  are  assumed. 

4.1.  Logarithmic  Quantizers 

A  logarithmic  quantizer  is  modeled  as  a  base-2  logarithmic 
nonlinearity,  followed  by  a  uniform  quantizer.  In  this  way,  an  M-level 

k 

quantizer  has  the  threshold  spacings  distributed  at  *  al  for  k 
representing  integers  between  -  (m/2-1)  and  (M/2-1).  In  this  case,  or  is 
denoted  as  the  threshold  spacing  for  the  logarithmic  quantizer.  Note 
that  this  definition  is  different  from  the  one  defined  by  Tong  [15]  for  the 
uniform  quantizer.  For  the  computation  of  Pq(t)  in  this  case,  kor  in 

Eq .  (3.34)  is  replaced  by  ar2  .  Hence  p  (2 )  for  the  logarihmic  quantizer  is 

Q 

given  by 


y2> 


TTJ. 


E 

n=l 
n  odd 


M/2-1 
k=0  * 


ex 


p{- 


(2ka)2 


*n-l(2«) 


P  n(2) 
x 


(4.1) 


where  again  A^  *  \  for  k  =  0,  A^  =  1  for  k  #  0. 

4.2.  Analysis  of  the  MTI  Performance  with  Logarithmic  Quantizers 
i)  Single  Canceller  MTI 

Apply  Pq(2)  of  Eq.  (4.1)  to  y  of  Eq.  (3.32),  the  loss,  incorporated 
with  a  16  level  quantizer,  as  a  function  of  for  different  values  of 
a  is  plotted  in  Fig.  4.1.  The  values  of  p^  shown  are  ranging  from  0.1  to 
0.999.  As  p^  approaches  unity,  the  loss  is  infinite.  This  is  the  case 
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when  x(t)  =  x(t+D)  for  perfectly  correlated  input  clutter  ,  corresponding 
to  a  zero  MTI  output.  It  is  seen  that,  for  this  case,  as  the  input  clutter 
becomes  more  highly  correlated,  the  loss  rises  significantly.  The  quantizer 

input  correlation  function  p  (t)  is  thus  smaller  than  the  linear  input 

* 

correlation  function  PX(T)*  However,  all  these  are  justified  by  the  fact 

that  the  quantizer  tends  to  decorrelate  the  clutter  input  because  of 

the  spread  out  of  the  clutter  spectra  in  the  presence  of  a  quantizer. 

Hence  the  effectiveness  of  an  MTI  is  reduced  drastically.  In  case  of  lowly 

correlated  input  clutter,  for  instance  <  .9,  if  the  threshold  spacing 

is  not  too  big  or  too  small,  the  performance  of  the  digital  MTI  and 

linear  MTI  are  almost  the  same  regarding  the  degradation  in  performance. 

As  shown  in  Fig.  4.1.,  the  degradation  is  less  than  2.5  dB  for  p  <  .9. 

x 

As  a  result,  we  will  be  concentrating  more  on  the  performance  of  an  MTI 
for  highly  correlated  input  clutter,  in  the  sense  that  p^  2  0.9.  It  can 
be  seen  that  the  choice  of  at  has  a  dramatic  impact  on  the  resulting  loss 
of  the  digital  MTI.  Too  large  or  too  small  a  choice  of  at  (i.e.  or-**  or  ar-»0) 
will  give  rise  to  a  maximum  loss  of  4.58,  9.55  and  14.54  in  dB  for  input 
correlated  coefficient  =  .9,  .99  and  .999  respectively.  Consider  the 
8  and  32  level  logarithmic  quantizers.  The  loss  as  a  function  of  a >  is 
plotted  in  Figs.  4.2  and  4.3.  For  any  value  of  M,  it  is  found  that  there 
is  always  an  optimum  spacing  for  each  input  correlation.  Nevertheless, 
at  the  setting  of  optimum  spacing  for  both  quantizers  taken  for  p^  =  .999, 
the  minimum  loss  are  7.03  and  5.97  in  dB  for  the  8  and  32  level  quantizers 
respectively.  Hence,  the  difference  is  only  1.06  dB  in  terms  of  improvement. 
With  more  quantizer  levels,  the  loss  can  be  diminished  by  a  large  amount. 

For  instance,  when  p^  =  .999,  an  improvement  of  5  dB  is  achieved  for  a 


Figure  4.2.  Loss  of  single  canceller  MTI  with  8-level  logarithmic  quantizer. 


32 -level  quantizer  as  compared  to  a  8- level  quantizer,  even  though  the 

.9 

spacing  a  =*  2  is  nonoptimum.  Table  4.1  shows  the  minimum  loss  and 
optimum  spacing  with  different  number  of  output  levels  for  ox  *  *9,  .99 
and  .999.  It  is  found  that  a  16-level  quantizer  is  almost  as  good  as  a 
32 -level  quantizer  at  optimum  setting  of  threshold  spacing.  Although  the 
optimum  spacing  for  M  =  32  and  64  are  the  same,  it  can  be  shown  that  they 
vary  a  great  deal  at  setting  of  nonoptimum  spacing, 
ii)  Double  Canceller  MTI 

Again  the  loss  as  a  function  of  p^  is  plotted  in  Fig.  4.4.  Here  the 
loss  is  much  greater  than  the  loss  of  a  single  canceller.  Theoretically, 
there  is  a  progressive  translation  of  power  toward  higher  frequencies, 
and  therefore  the  high  filtering  effect  of  the  cancellation  circuits  is 
accentuated.  This  explains  why  the  loss  for  the  double  canceller  is  so 
much  higher  than  that  for  the  single  canceller.  The  loss  for  8-  and 
32-level  logarithmic  quantizers  is  obtained  and  plotted  in  Figs.  4.5 
and  4.6.  The  result  shows  that  for  spacing  less  than  optimum  spacing,  the 
loss  will  be  differentiated  by  a  wide  margin  for  M  *  8  and  32;  while  for 
spacing  greater  than  the  optimum  spacing,  the  loss  is  within  1  dB. 
Nevertheless,  in  the  limit  when  the  quantizer  spacing  becomes  too  large 
or  too  small,  the  loss  can  1 »ach  10.32,  24.83  and  39.78  in  dB  for  =  .9, 
.99  and  .999  respectively.  Table  4.2  shows  the  minimum  loss  and  optimum 
or  for  different  number  of  output  levels. 

4.3.  u,-Law  Quantizers 

A  u-law  quantizer  is  represented  by  a  memoryless  nonlinearity  function 


g(x)  =  V  log(l+ux/V) /log(l+u)  (where  u  >  0  and  V  >0),  followed  by  a 
uniform  quantizer.  The  nonlinearity  is  defined  as  a  u-law,  or  logarithmic 


Figure  4.4.  Loss  of  double  canceller  MTI  with  16-level  logarithmic  quantizer. 


compression,  curve  [24]  where  |i  is  a  compression  parameter  which  controls 

the  degree  of  compression  and  may  be  chosen  so  that  large  changes  in  the 

input  produce  relatively  small  charge  in  the  output;  and  V  is  the  preset 

maximum  input  voltage.  The  nonlinearity  spreads  out  low  amplitude  input 

clutter  over  a  large  range  while  shrinking  the  higher  amplitude  input 

clutter  into  a  smaller  region.  The  effect  is  to  allocate  more  quantization 

levels  to  the  lower  amplitudes,  which  generally  have  higher  probability. 

In  this  case,  the  spacing  of  the  quantizer  levels,  which  is  a  function  of 

V  k/V  k 

u,  and  V,  is  located  at  a.  -  -  [(1+pi)  -11  instead  of  or2  as  for  the 

K  Pi 

logarithmic  quantizer.  Apply  a^  to  Eq.  (3.34),  Pq(T)  f°r  the  H-law 
quantizer  is  given  by 


,  <T)  -  Z  f  M/|'V  eKp[-£t<l+u)k/’'-U12/2lH  [2[(nu)k/v-l]) 

A  TTOn  n=l  L  k=0 
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where  4  ■  j  for  k  *  0,  A,  =1  for  k  5*  0. 

K  tC 

4.4.  Analysis  of  the  MTI  Performance  with  u-Law  Quantizers. 


i)  Single  Canceller  MTI 

With  p.  ■  1,  10  and  100,  V  =  100,  the  loss  is  plotted  with  respect  to 
from  .1  to  .999  for  a  16 -level  quantizer  and  is  shown  in  Fig.  4.7. 

As  analogous  to  the  logarithmic  quantizer,  the  loss  of  the  digital  MTI 
increases  significantly  as  the  input  clutter  becomes  highly  correlated, 
i.e.  for  ox  2  -9.  The  quantizer  experiences  the  biggest  loss  when  a 
compression  parameter  of  100,  the  largest  of  the  three,  is  chosen  for 
^  .9.  When  u.  is  large,  the  output  levels  are  crowded  about  zero. 
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In  the  limit  when  p.  approaches  +«  ,  the  quantizer  behaves  essentially  as 
a  hard  limiter,  the  maximum  loss  is  exactly  the  same  as  that  obtained 
for  the  uniform  and  logarithmic  quantizers.  For  highly  correlated  input 
clutter,  u  =*  100  comes  out  to  be  the  best  of  the  three.  Figure  4.7  does 
indicate  that  the  choice  of  p.  is  vital  to  the  performance  of  digital  MTI. 
Again,  the  following  discussion  will  be  restricted  to  highly  correlated 
input  clutter  with  £  0.9.  Results  of  8  and  32  level  quantizers  with 
V  =  10  are  shown  in  Fig.  4.8  and  Fig.  4.9.  Intutively,  for  any  given  M, 
there  is  an  optimum  u  which  gives  rise  to  a  minimum  loss.  Indeed,  the 
optimum  u  occurs  for  each  p  .  However,  it  is  interesting  to  find  that  the 
optimum  u  tends  to  increase  as  the  input  correlation  function  increases 
As  optimum  p's  of  M  =  8  are  compared  to  those  of  M  =  32,  it  can  be  seen 
that  with  more  quantizer  levels,  the  optimum  compression  parameter  is 
incremented  by  a  great  deal.  For  instance  at  p^  =  .999,  the  compression 
parameters  at  minimum  loss  are  79  and  165  for  M  *  8  and  32  respectively. 
With  V  set  at  1000,  the  loss  of  8  and  32  level  quantizers  are  plotted 
in  Figs.  4.10  and  4.11.  The  results  here  are  analogous  to  the  V  *  10  case 
Again  each  of  the  loss  curves  shows  an  optimum  j,  coupling  to  a  minimum 
loss.  A  close  examination  of  Figs.  4.8-4.11  reveal  some  interesting 
results.  Each  of  the  optimum  p.'s  and  their  corresponding  loss  at  V  =  1000 
is  smaller  than  that  obtained  at  V  *  10.  It  implies  that  an  improvement 
in  the  performance  can  be  achieved  by  a  combined  action:  increasing  the 
number  of  quantizer  levels  and  allowing  a  large  preset  input  voltage. 

Table  4.3  lists  the  optimum  y,'s  and  their  corresponding  loss  with  V  =  10 
for  various  output  levels. 


Compression  Parameter  p 

igure  4.8.  Loss  of  single  canceller  MTI  with  8-level  U'law  quantizer  and  V  =  10. 


Compression  Parameter  pi 

Loss  of  single  canceller  MTI  with  32-level  p,-law  quantizer  and  V  =  1000. 
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ii)  Double  Canceller  MTI 

Figure  4.12  shows  the  loss  as  related  to  p  which  goes  from  0.1  to 
0.999.  The  loss  curves  here  produce  the  same  trend  as  those  obtained  for 
the  single  canceller  MTI,  but  the  loss  grows  much  higher  than  the 
single  canceller  case.  Figures  4.13  and  4.14  show  the  loss  as  a  function 
of  p  with  8  and  32  level  quantizers  for  the  case  V  =  10.  It  is  seen 
that  the  loss  is  insensitive  to  the  value  of  u  chosen  since  each  curve 
is  rather  flat  as  compared  to  that  of  the  single  canceller  case.  When 
V  is  increased  to  1000,  the  loss  curves  are  done  over  and  the  data  is 
plotted  in  Figs.  4.15  and  4.16.  Here  the  loss  becomes  more  sensitive  to 
the  value  of  p.  For  all  the  figures  shown,  an  optimum  a^  can  always  be 
found.  Finally,  more  quantizer  levels  and  larger  preset  input  voltage 
are  favorable  towards  the  optimum  performance  of  a  double  canceller  MTI. 
Again,  there  is  always  a  limit  on  the  loss  here.  Table  3.4  lists  the 
optimum  p's  and  their  loss  with  V  =  10  for  values  of  M  ranging  from  4  to  64. 
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Figure  4.12.  Loss  of  double  canceller  MTI  with  16-level  n-law  quantizer 


Figure  4.15.  Loss  of  double 
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5 .  CONCUJSIONS 

It  has  been  shown  that  quantization  effects  on  digital  moving  target 
indication  depend  on  the  number  of  output  quantization  levels  and  the 
spacing  of  these  levels.  Among  the  quantizers  discussed  before,  the 
clutter  attenuation  is  being  used  as  the  performance  criterion.  It  has 
been  shown  that,  in  general,  uniform,  logarithmic  and  p.-law  quantizers 
give  almost  identical  performance  for  a  few  quantizer  levels.  As  a 
matter  of  fact,  a  4-level  uniform  quantizer  is  equivalent  to  a  4-level 
logarithmic  quantizer.  As  the  number  of  quantizer  levels  increases,  the 
uniform  quantizer  gives  the  best  performance  and  logarithmic  quantizer 
presents  the  worst.  However,  as  the  comparison  parameter  approaches 
zero,  a  p.-law  quantizer  will  be  similar  to  a  uniform  quantizer.  Moreover, 
it  is  likely  to  happen  that,  with  a  large  enough  number  of  output  levels, 
the  loss  will  not  be  improved  by  any  significant  amount.  Hence,  it  is 
concluded  that  the  uniform  quantizer  is  the  best  fitted  in  the  sense  that 
the  performance  of  digital  MTI  using  the  uniform  quantizer  gives  a 
minimum  loss  corresponding  to  an  optimum  spacing,  under  the  assumptions 
that  the  input  clutter  process  is  Gaussian  and  the  quantizer  has  uniform 
step  size.  Finally,  it  has  been  verified  that  there  exists  a  common 
bound  on  the  loss  for  uniform,  logarithmic  and  p.-law  quantizers. 

It  was  seen  that,  in  general,  uniform,  logarithmic  and  p.-law 
quantizers  give  almost  identical  performance  at  few  quantizer  levels. 

As  a  matter  of  fact,  a  4-level  uniform  quantizer  is  equivalent  to  a 
4-level  logarithmic  quantizer.  As  the  number  of  quantizer  levels 
increases,  uniform  quantizer  gives  the  best  performance  and  logarithmic 


quantizer  presents  the  worst.  However,  as  the  compression  parameter 
approaches  zero,  a  ti-law  quantizer  will  be  similar  to  a  uniform 
quantizer.  Moreover,  it  is  likely  to  happen  that,  with  a  large  enough 
number  of  output  levels,  the  loss  will  not  be  improved  by  any 
significant  amount. 
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